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THE  SQUARE  ROC?  X2TH0D 

Prepared  by t 
Calvin  C.  El got 


ABSTRACT:  Condition®  under  which  solution  of  a set  of  linear  equations  can  be 

effected  by  the  square  root  aethod  are  given.  While  noch  literature  exists  on 
the  technique  of  this  aethod,  we  have  bean  unable  to  find  a discussion  of  con- 
ditions under  which  this  technique  succeeds.  A aethod  of  decomposing  a large 
class  of  hendtion  aatrioes,  (soae times  requiring  quaternions),  also  appears  to 
be  new. 
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IT  October  lff2 


This  report  contains  the  result  of  an  investigation  carried  out  la  connection 
with  a problem  in  least  square  fitting  of  a function  which  expresses  ocaplex 
jaw  ef  a spin-stabilised  projectile  la  tens  of  distance  down  range*  This 
problaa  was  resolved  Into  a system  of  linear  equations  with  complex  coef- 
ficients, the  solution  of  which  is  discussed  here*  The  work  was  carried  oat 
under  the  sponsorship  of  (MR  project  IR-Q44-003»  Maaerlcal  Analysis  and  BuOrd 
project  X0L-Se3d-452-l-53,  Tree-Flight  Aerobellistios  of  Spin-Stabilised  Bodies, 


EDUARD  L.  WOODKARD 
Captain,  USB 
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THE  SQUARE  ROOT  METHOD 


1.  Introduction 

Throughout  this  paper  S * (s.  j)  will  denote  a triangular  matrix, 

(i.e.,  1 > J implies  s^j  « 0),  *A  will  denote  the  transpose  of  the  matrix 

>*:  R,  C,  Q will  denote  respectively  the  field  of  real  numbers,  the  field  of 
complex  members  and  the  afield  of  quaternions  with  real  coefficients.  M(F) 
will  denote  the  ring  of  nxn  matrices  with  elements  in  a sfield  F. 

In  the  square  root  sethod  one  solves  the  system  of  wjnatiwna 

AX  • G for  X, 

(A  * MjjCr),  X en  unknown  column  vector,  Q a given  column  vector),  in  the 
following  merger  (of.  ref.  (a))t 

One  finds  s matrix  S such  that  *S.  S * A,  Assuming  A non- singular,  it 
t 

follows  S and  S are  non-singular  so  that 

VSX  « 0 implies  SX  « X 

t -1  t 

where  I * S .0.  Because  of  the  simple  forms  of  S and  £,  the  calculation 

of  %"*.  G end  of  X ie  fairly  simple.  For  another  diaomssien  of  this 
mothod  and  further  referenosa  see  reference  (b). 

2.  The  Theorems 

In  what  follows  ws  are  given  a sfield  F and  an  anti-automorphiam/of  the 

afield  such  that  1.  If  B ■ (b^j)  * Kj,(F)  and  B * (<rb^j)  then 

6 *B  *Vb.  He  shall  sometimes  interpret  B * M_(F)  as  s linear  trsnsfoi mstion 
on  the  left  n-dimcdslonal  vector  spaoe  of  n»tup£es  of  F over  F.  To  say  B is 
non-singular  means  B has  on  inverse  which  is  equivalent  to  the  requirement  that 
the  1st  >•  of  B be  the  whole  spaoe.  This  latter  condition  is  easily  seen  to 
equivol  r.  to  the  requirement  that  the  rows  of  the  matrix  B be  left  linearly 
independent. 

i-*"”—  a triangular  matrix  S « M^F)  is  non-singular  if  and  only  if 
=ii  ^ °»  *cr  * * *»  * ’ » »• 

(1)  For  n * 1,  the  lemma  is  trivial.  Suppose  Tf  / 0,  and  n > i. 

We  nay  assume  inductively  that  the  matrix 
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[ ll 

*22  8 23  •••  *2n 

0 ®33  **•  *3n 

• It  • • • • 

0 0 ...  .0  a^ 

has  left  linearly  independent  rows.  It  the  row  vectors  of  8 ere  dependent  then 
there  exists  b « F ouch  that  b e^j  * 0.  Since  sjj^  0,  b « 0*  Thie  implies 

the  dependence  of  the  row#  of  the  netrix  £l]  which  contradicts  the  inductive 
assumption. 

low  suppose  TT  * 0.  If  By  ■ 0 for  sons  1 > 1 then  hy  induction 
the  rows  of  f.x"}  we  dependent  end  hence  the  z-ovs  -f  S«  I-  j|  "£s_  J*  than 
a^  ■ 0.  The  n vectors  (s^,  *13  • • • ■ic)  * 1 - 1 - *»  dependent  since 

they  lie  in  an  r-1  diwmsional  spaee.  If  2!  t*i^*i2>  •••  * *in^  * ® then 

/v  / 

^i'*!!*  ®i2»  8 in'  " 

/ »/ 

If  A * M^F),  then  for  1 £ k S n,  A*  c %(F)  denotes  the  natrix  obtained 
fro*  A by  delating  all  rows  and  eolunna  except  rows  and  ooltmms  1 through  k. 

Theorem  1 

Suppose  F commutative.  Let  A * 74  (F)  be  non- singular.  If  there  exiete 
S 0 M^F)  such  that  <*"**5.  S * A then 

(a)  /Aj/  ^ 0,  1 * 1,  — , n;*/^i/"  denotes  the  determinant  of  Ai. 

(b)  <rx, . x,  » /%/»  i * 1,  — , n has  a eolution  in  F. 

(c)  <^A  A. 

£xxja 

Sinoe  /**$/•  /S/  * /A/  and  A is  non-singular,  /S/  0 and  //rS/  0. 

Simn  /S/  * TT  aljL,  s^  / 0 and  hence  ^(b^)  T 0.  It  followa  that 
/A^/'  0,  =’l,  — , n - 1,  since  <f  *3^  « A^,.  Thus  (a)  ia  satisfied. 

For  n - k 1,  we  have 

0 f /^J  + /4k-i/  “ *■  *ii-*u/  TT  <r®U**ii  “ /ftkk*ekk 
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For  k ■ 1 we  have  6 Sj^.  s^  * a^.  Ase^cdsf  fx^.x^  » /Aj^/  w ***** 

<r*o'3V  ^•kk*  *kk  * ^*0  *kk )'  • x©  rj^  ■ /Ajp/i  (b)  la  demonstrated. 
To  demonstrate  (o)t 


V /l 


V C <f,r 


.si  - rCs.*ff)-  v* 


Jhferaa  2 

lot  F be  a afield.  Let  G denote  the  aet  of  x e T such  that  <x  * x* 
Suppose  for  every  g ( G there  exists  x < F such  that  fx,  x * g.  Giron 
a poa-alngular  natrix  A « M^F).  If 

<a)  A^,  i * 1,  ...  , m is  non-alngalar. 

(b)  V A * A (i.e.,  / a^  ■ a^). 


Than  there  exists  8 « IL^F)  such  that  V S.S  - A. 


Proof 


We  observe  first  the  eedatanoe  of  auoh  an  8 is  equivalent  to  the  existence 
ef  a1  j *F,  l-i-J-n»  auoh  that 


[2]  < 

t '•rl-  *rj  * *ij*  J 21  1 ,ine* 

f A'***-'  **J  ’ </  srJ-  *rl  * ^*11  " »Ji- 

For  n * 1 the  theoren  la  trivial*  We  note  that  conditions  (a)  and  (b)  hold 
for  (with  (n— l)  replacing  a)*  We  asraaa  Inductively  the  existence  of 

a1j,  1 i i ^ J < a auoh  that  [2]  is  satisfied.  We  seek  now  a^  i « 1,  ...» 
satisfying 


$■  f a ..  a_  “ a . 
ri  rn  in 

A»t 


n 


For  i < n,  we  need  solve  a r equation  the  coefficient  of  the  unknown  of 
which  is  bh-  Since  A^  la  non-singular  the  leas*  tells  ua  a^  ^ 0 and  so 

o' en  / C Xt  remains  only  to  find  iM  metis fyiag 


£ <T 

si  ■=/ 


* a_ 


We  observe  that  <r  a^  • *^5  m ( 


’rn 


rn 


’rn 


) ■ o'  a, 


m 


rn 
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Therefore  ( 6 *n  . s^)  G.  Since  - au,  /a^  so  that 

tjjjj  c G.  Cue  sees  easily  that  0 is  ©a  abelian  group  ©o  that 

Sm  - ^ ^ru  • 8rn  ie  ln  0 

/»■/ 

and  the  result  follows. 

3.  The  Special  Casas 

Corollary  1.  L©t  A ( M^C)  be  non- singular  and  syaaetrlo.  A necessary  and 
sufficient  condition  that  there  exist  S £ >^(0)  such  that  *S.  8 * A is  that 
/Aj/,  i * 1,  •.  , n - 1,  be  non-sero. 

Let  F * C and  ^ « 1. 

Corollorr  2.  Let  A £ M^(Q)  be  non-singular.  Suppose  Va  * A where 

<T  (*!  + *2J)  * ix  + s2J,  (^  - a + bi,  »2  - o ♦ di). 

A nocessary  and  sufficient  condition  that  there  exist  S £ M^Q)  such  that 
S.S  * A is  that  A^,  i *1,  ...  , n - 1 bo  non-singular. 
fiCSfil 

One  sees  readily  6 is  an  anti-autoaorphisa  of  Q.  G consists  of  the 
eleaents  of  the  for*  b + s j with  b real.  For  b,  s fixed  we  seek  s^  and 

such  that  (s^  + *2J)(*i  ♦ *2^)  ■ b + * J.  If  s * 0,  this  is  easy.  Assua© 
t fl  0.  s^,  Sj  auet  be  non-sero  and  aust  satisfy  s^ij,  - *2*2  * b and  2s^s2  » s. 

This  is  equi Talent  to  s2  ■ 7—  and  s^i^  - (si/As^s^)  * b.  The  latter  is 

equivalent  to  A(*j*i)^  - Abs^i^-si  * 0 which  always  has  a solution, 

via.  ■ 2 (h  + \jb2  + si)  > 0. 

Ve  euppleaent  oorollory  2 with  the  renark  that  if  A £ M^(Q)  is  heraitien, 

and  if  or  chooses  the  diagonal  eleaents  of  S properly,  (identifying  * t D J with 
*),  S co!  ains  eleaents  of  a particularly  slaple  fora.  If  we  call  s + 0 J a 
pon»l«  sm&nla*  end  0 <«■  s j a ***»<*  that  every  eleaent  of  8 ie 

either  ooaplex  or  chaste.  Indeed  the  Sj^  nay  always  be  chosen  still  aore  slxply 
as  either  a positive  real  nuaber  or  a positive  real  maber  tines  j.  These 
reaarks  aay  be  verified  using  ths  equations 


A 
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< — > 

*ii  a (aU  " 2 ^ "ri  * **t  i > 1 

A~t 

Bij  * *aiJ  “ 2<r*rl  * Brp/*li  ^ > 1 >2^ 

• A*  i 

«n  " (‘u)  > ‘is  * *iy'au 

and  tin  facta  s the  inverse  of  a ohk.ste  *oataraA«i  la  ohaate,  the  product  of  a 
and  chaste  quaternion  la  chanie,  crg^  . a^  1*  ooaplex  whether  a^  !• 
ooarploDC  *r  abasia* 

Gorollory  2 and  3 insure  that  a poaltlTe  definite  re*l 
a poeitive  dafinlta  herdtian  Matrix  are  deooapooatl®  in  the  esapisx  asd  ^uaierdon 
Wm  respectively.  However  in  these  oaaaa  tha  daaoapoaltioa  la  possible  in  tha 
rael  oeer.pl  *a  d&aal&s  respectively*  ?h*.a  la  basause  a^j*  * (positive  reel  wssber) 

ass  « salntlrta  in  the  reals  as  does  f m*t  » «u  ■ (positive  real  unbar).  For 
be  interesting  discussion  of  herdtlan  uatrioaa  aaa  rafaranoa  3. 

Utilising  theorem  2 and  "Tha  Principle  of  tha  Irrelevance  of  Algebraic 
laeqaali ties , * ( raf.  4#  P*  4),  one  oan  prove  that  tha  naoaasary  conditions 
asnaciated  in  theoran  1 for  tha  axlatanoa  of  a daoaapoaitlon  of  a non-aingalar 
mtrix,  (in  the  case  of  a coaamtative  F) , la  Indeed  suffloiar/c.  Tha  raenlta  of 
tha  above  paragraph  ara  than  obtainable  directly  from  this  theorssu 
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